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ABSTRACT 


The  errors  are  anaJ-yzed  which  arise  when 
using  ICasemir's  method  to  measure  air-earth 
current -density.  Step  function,  ramp  function, 
and  periodic  functions  in  the  current  as  well 
as  in  the  time-constant  theta  of  the  air  are 
discussed. 
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THE  EFFECT  OF  MISMATCHING  IN  THE  MEASUREMENT 
OF  THE  AIR-EARTH  CUEREIW-DEI^SITY 


INTRODUCTION 

Successful  measurements  of  the  air-earth  current  density  started  in  I906 
when  Wilson!  introduced  his  well-known  method  which  was  based  on  time  sampling 
and  was  not  directly  suitable  for  continuous  recording.  The  difficulty  in 
making  a  continuous  recording  waa  the  effect  of  the  potential  gradient  on  the 
current  measurement.  In  1933  Scrane^  tried  a  method  of  compensating  for  this 
effect,  but  he  did  not  really  succeed  although  his  method  is  feasible.  Adam- 
son3  rearranged  the  method  of  Scrase  and  used  a  modem  field  ttiiII  to  over¬ 
come  the  sluggishness  of  the  polonium  collector  used  by  Scrase.  A  much  simp¬ 
ler  circuit  for  continuous  c\irrent  measurement  was  introduced  in  1955  ty 
Kasemir,^  who  conqpensated  for  the  displacement  current  generated  by  fluctua¬ 
tions  of  the  field  in  an  ingenious  way.  The  errors  of  this  method  and  the 
limitations  for  practical  use  are  discussed  in  detail  in  this  report. 

The  three  electric  parameters  most  often  used  to  determine  atmospheric 
electric  occurrence  in  fair  weather  are  the  field  strength  F,  the  sum  of  the 
polar  conductivities  A  ,  end  the  conduction  current  density  i.  They  are  com¬ 
bined  in  the  Ohm's  law  which  has,  in  the  three-dimensional  case,  the  fora 

i  =  A  •  F.  (1) 

Normally  it  is  sufficient  to  continuously  record  two  of  these  quantities  in 
order  to  determine  all  three.  However,  since  the  measiu:ing  instruments,  and 
sometimes  the  method  used,  give  the  measxired  parameters  with  a  more  or  less 
significant  error,  it  shoxild  be  noted  that  a  measurement  of  three  vari¬ 
ables  is  necessary  to  have  a  continuous  check  of  proper  functioning. 

One  of  the  most  significant  errors  comes  from  the  influence  of  the  dis¬ 
placement  current  in  air  on  a  conduction -ctirrent  measurement.  This  report 
shows  theoretically  the  expected  error  when  using  Kasemir's  method  of  compen¬ 
sation.  Starting  from  the  general  dif f erentiaG.  equation,  two  cases  will  be 
treated  separately.  One  is  the  case  when  all  fluctuations  originate  from 
changes  in  the  current  density  and  the  other  is  the  case  when  nil  fluctuations 
come  from  changes  in  the  conductivities.  It  is  possible  to  treat  these  vari¬ 
ables  separately  because  the  current  and  the  conductivity  can  be  taken  in 
approximation  as  independent  variables. 

In  discussing  the  solution  of  differential  equations,  it  is  practical  to 
consider  given  time  functions  for  the  variables.  A  step  function  and  a  ramp 
function  are  chosen  because  they  demonstrate  the  mismatching  condition  ade- 
q\iately.  A  step  function  is  chosen  because  it  represents  sudden  changes  in 
the  variables;  and  a  ramp  function,  which  is  a  function  increasing  llnearily 
with  time,  is  chosen  to  represent  slow  changes .  In  addition,  a  frequency- 
amplitude  discussion  for  periodic  fluctiiations  of  the  variables  is  set  forth. 
Frequency-anqjlitude  diagrams  in  a  double  logarithmic  scale  as  introduced  by 
Bode5  are  used  to  interpret  the  results. 


DISCUSSION 

Derivation  of  the  Differential  Eqmtion 

Kjasemir's  metliod  of  current  measurement  consists  of  an  insulated  antenna 
that  is  exposed  to  the  atmosphere.  This  current  sensor  is  grounded  by  a 
resistor  with  a  suitable  capacitor  in  pai:allel.  The  current  through  the 
resistor  is  measured  by  a  sensitive  current  meter.  The  necessary  conditions 
for  a  true  measurement  of  the  conduction  current  are:  the  measuring  resistor 
is  so  small  that  the  voltage  drop  over  it  is  snail  compared  with  the  open 
loop  voltage  at  the  antenna,  and  the  time  constant  T  of  the  EC  circuit  is  at 
1  times  equal  to  the  time  constant  of  the  air.  In  practice  these  conditions 
are  fulfilled  in  approximation  only,  even  for  measurements  at  the  ground. 
Special  difficulties  arise  when  current  densities  are  to  be  measured  as  a 
function  of  height  with  balloon  soundings  because  the  conductivity  of  the  air 
is  changing  exponentially  with  height.  Kasemir,^  who  made  the  first  current 
measurements  with  balloons,  has  EB.de  a  detailed  analysis  of  the  measuring 
technique  involved. 

For  the  antenna,  the  equivalent  circuit  diagram  as  given  by  Kaseaiir  and 
Ruhnke”^  is  used.  It  is  shown  in  Fig.  1  where  F  is  the  field  strength,  h  the 
effective  height  of  the  antenna,  c  the  antenna  capacity,  and  r  the  internal 
resistance  of  the  antenna.  The  product  rc  is  equal  to  the  time  constant  9  of 
the  air  surrounding  the  antenna.  The  condition  that  the  measuring  resistance 
R  should  be  small  is  equivalent  to  the  condition  that  the  current  measured  is 
the  short  circuit  current  in  the  diagram  of  Fig.  1.  This  is  normally  ful¬ 
filled  to  a  satisfactory  extent  with  a  measuring  resistance  of  10^  ohm  because 
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Fig.  1.  The  equivalent  circuit  diagram  for  measuring  current  density. 
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the  internal  resistance  r  of  the  current  antenna  will  he  in  the  range  of  10^3 
ohms  when  measured  at  the  gr  ind.  A  more  elegant  way  to  satisfy  this  condi¬ 
tion  lies  in  the  use  of  the  analog  computer  technique.  An  operational  ampli¬ 
fier  TTith  an  electrometer  input  with  the  SC  network  in  the  feedback  loop 
replaces  the  current  meter  and  the  SC  network  of  Kasemir's  method.  This 
arrangement  is  especially  practical  if  current  densities  are  to  be  measured 
in  balloon  and  rocket  soundings  since,  with  the  ejqponential  increase  of  con¬ 
ductivity  in  high  altitudes,  the  short  circuit  condition  is  easier  fulfilled 
than  -^-rith  the  traditional  arrangement.  The  equivalent  circuit  diagram  for 
this  method  is  shown  in  Fig-  2  where  the  gain  A  of  the  operational  amplifier 
for  convenience  is  assumed  infinite.  Practical  values  of  A  are  between  60  db 
and  120  db.  It  can  be  shown  that  the  effective  input  impedance  of  this  cir¬ 
cuit  at  the  input  of  the  amplifier  is  1/A  times  the  feedback  impedance. 
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Fig.  2.  Circuit  with  operational  amplifier  to  measiure  current- density. 

One  more  definition  is  necessary  to  set  up  the  general  differential  equa¬ 
tion  for  this  circuit;  namely,  the  desired  conduction  current  density  i  and 
a  factor  'A  which  has  the  dimension  of  an  area  and  means  the  effective  area  of 
the  antenna.  The  open-loop  voltage  Fh  is  related  to  the  conduction  component 
Mi  of  the  short-circuit  current  by  the  internal  resistance  r. 

Fh/!iL  =  r.  (2) 

The  differential  equation  governing  this  circuit  is  then  given  by 

U/R  +  C  dU/dt  =  -(Fh/r  +  c  d?h/dt).  (j) 

With  equation  (2)  and  with  the  relation  for  the  time  constant,  the  differ¬ 
ential  equation  can  be  rearranged. 
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U/R  +  C  dU/dt  *  -(I4i  +  M  di9/dt)  . 


(4) 


Tlois  eq\iation  is  valid  for  homogeneous  conductivity  and  current  density  in 
the  neighborhood  of  the  entenna.  The  current  as  well  as  the  conductivity  w&y, 
however,  be  functions  of  time.  The  general  solution  with  no  initial  charges 
assumed  may  be  expressed  by 

t 

U  =  -Mir  -  exp(-i)  I  i  ^  +  (9  -  T)  ^  exp  (i)  dt.  (5) 

T  T  o  _  dt  dtj  T 

It  can  be  seen  that  for  9  =  T  =  constant,  the  so-called  matching  condition, 

U  is  proportional  to  the  current  density  i  for  any  time  function  of  i.  This 
is  no  longer  the  case  if  the  conductivity  changes.  Then,  d9/dt  is  no  longer 
zero  and  this  term  introduces  an  r  error  term  into  the  measurement.  But  in 
addition  9  is  no  longer  equal  to  T,  and  any  current  change  di/dt  brings  an 
additional  deviation  from  the  time  measurement.  The  (l)  solution  of  equation 
(4)  for  the  time  functions  in  question  will  be  obtained  by  using  the  Laplace 
trans  f o  rmat ion . 

The  Step  Function  and  the  Ramp  Function 

To  demonstrate  the  influence  of  mismatching  and  conductivity  fluctuations, 
a  step  function  for  the  current  i  or  the  time  constant  9,  respectively,  is 
assumed.  In  the  first  case  let  9  =  constant  and  the  current  jump  from  the 
value  i  to  i  +  Ai  at  t  equal  zero.  The  Laplace  transformation  for  this  case 
is  given  by 

IlXii  +  sU(s)  C  -  ^  Ut=0  =  -Mi(s)  -  sM9i(s)  +  Mi^_Q  ©  (6) 

R  R  ”  ~ 


with 


l(s)  .  iUJi;  .  i; 


As  a  matter  of  convenience  let  the  equation  be  normalized  to  a  dimension¬ 
less  one  with  the  help  of  the  relation 


-MiR  ■ 


Low  the  Laplace  equation  can  be  easily  solved  for  V(s) 


V(s)  = 


U 


+  Ai  (9  -  T) 

1  (1  +  sT) 


The  time  relation  V(t)  will  be  obtained  after  bach  trans foruau ion  into  the 
time  domain.  It  contains  all  transients  for  uhe  given  initial  conditions  and 
represents  the  normalized  output. 

V(t)  =  1  +  I±:lL  exp(-t/T)  t  ^  0 

i  i  T 

V(t)  =  1  t-h  0 


t-h  0 
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(5) 
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This  e:':pression  means  -eneraJJLy  that  changes  of  the  current  which  occur  fast 
ccmpared  to  the  tirie  constant  T  are  either  exaggerated  or  damped  in  the 
recordea  output ^  depending  on  the  natching  condition. 

In  "ig.  3,  relation  (9)  is  plotted  against  tine  for  three  different 
matching  conditions;  9  =  T;  9  =  l/ZT;  9  =■  ‘27. 


Tircuit  response  zo  current -changes  for 
three  different  matching  conditions. 


irrors  cue  tc  changes  of  zhe  current  slow  to  the  tine  constant  T  can  he 
treated  r/  assuming  a  raitp  function  of  the  current.  Let  i  (t)  be  eq’ual  to 
+  di/ -it  t,  "ith  di/ it  the  rata  of  current  change  ter  time  unit.  The  Laplace 
ransfortation  then  is 

t  "  P  UO) 

-he  initial  conditions  are  given  by; 


=  i 

If  this  is  applied  to  equation  (6)^  and  V(s)  is  solved  for, 

V(s)  =  i  +  L.  _^  +  i  (Q-T)  ^  ^  T(T-9)  ^ 

s  g2  Ati  s  i  At  (1  +  sT)  iAt 

In  the  time  domain  this  equation  reads: 

V(t)  =  1  +  ^  i  ^  ^  exp(-t/T) .  (12) 

At  i  i  At  i  At 

The  first  two  terms  represent  the  true  current  component.  The  term 
(9-T)  Ai/At  is  an  error  component  and  means  that  the  current  measurement  is 
changed  in  an  amomt  proportional  to  (9-T)  and  to  the  rate  of  current  change. 
The  last  term  represents  the  transient.  Summarizing,  fast  current  changes  are 
exaggerated  in  the  i*eading  by  the  factor  (9-T)/T  and  slow  changes  bring  a  con¬ 
stant  error  of  the  amcunt  (9-T)  Ai/At. 

In  Fig.  k-  relation  (12)  is  plotted  to  show  the  influence  of  slow  current 
changes.  Again  the  following  matching  conditions  are  used:  9  =  T,  9  =  1/2T, 
and  9  =  2T. 


In  the  next  case  let  i  be  constant  and  the  time  constant  9  of  the  air 
change  at  t  =  0  in  a  sudden  jump  from  the  value  9  to  9  +  A9.  Equation  (5) 
can  then  in  its  normalized  form  be  expressed  for  this  special  time  function  by 

V  =  1  +  a9/T  exp(-t/T)  t  0  (13) 

V  =  1  t  0 

Values  of  T,  1/21,  and  21  are  assiimed  for  A9  and  the  cutout  is  plotted  in 
Fig.  5. 

It  is  interesting  to  notice  that  the  error  term  of  eqiiation  (I3)  is  no 
longer  zero  for  9  =  T,  but  completely  independent  from  the  case  of  matching. 

The  constant  part  of  the  time  constant  9  does  not  even  appear  in  equation  (13 )• 
The  magnitude  of  the  error  is  proportional  to  the  change  of  9  and  inversely 
proportional  to  the  time  constant  T  of  the  measuring  device.  Changes  of  9, 
sTtial  1  compared  to  the  time  constant  T,  will  not  affect  the  current  reading 
appreciably.  In  practice  the  conductivity  has  surong  fluctuations  if  meas¬ 
ured  with  a  Gerdien  tube  of  conventional  size,  however,  the  current  antenna 
has  an  appreciable  area  and  the  mean  value  of  the  conductivity  over  this  area 
is  responsible  for  the  error  term  in  equation  (13).  In  practice  this  mean 
value  of  the  conductivity  shows  no  rapid  fluctuations,  and  the  error  term  is 
neglectible  if  current  antennas  with  effective  axeas  over  10  m2  are  used. 

For  slow  changes  of  9,  again  a  ramp  function  is  assumed  of  the  form 
9(t)  =  9  +  A9/At  t.  3y  using  the  Laplace  transformation,  the  normalized  solu¬ 
tion  of  equation  (4)  is  obtained  as 


c 


v(t)  =  1  +  ^  (1  -  exp  -t/T).  (l4) 

At 

Again  there  is  no  matching  possible,  and  the  error  is  proportional  to  the 
rate  of  change  of  the  time  constant  A©/ At. 

Periodic  Functions 


The  discussion  of  a  step  function  and  a  ramp  function  provides  a  good 
demonstration  of  the  effects,  but  is  not  sufficient  to  describe  the  actual 
behavior,  since  random  fluctuations  may  occur  in  the  current  as  much  as  in 
the  conductivities  i-rith  no  similarities  to  step  or  ramp  functions.  A  more 
precise  description  of  the  phenomena  involved  is  given  by  a  con^jlex  frequency 
characteristic  of  the  system.  There  are  often  quasi  periodic  fluctuations 
observed  in  the  recordings  of  atmospheric  electric  paxameters,  and  vith  the 
knovn  frequency  characteristic  of  the  system  it  is  easy  to  estimate  the  exist¬ 
ing  error  for  a  specific  frequency.  The  sine  wave  is  probably  of  all  time 
functions  the  one  which  best  approximates  the  measured  fluctuations.  Besides 
this,  it  is  possible  to  present  any  time  function  by  its  spectrum  obtained  by 
Fourier  analysis.  The  spectrum  of  the  output  is  then  equal  to  the  spectrum  of 
the  input  times  the  frequency  characteristic  V(jco)  of  the  system.  The  error 
term  of  the  system  is  given  by  the  deviation  of  the  transfer  characteristic 
from  unity.  In  other  words,  a  correct  current  measurement  is  obtained  if 
V(j<o)  is  equal  to  unity  for  all  frequencies. 

To  obtain  the  transfer  function  of  the  system,,  assume  periodic  functions 
of  i  and  9  and  consider  the  solution  of  equation  (4)  after  decay  of  the  transi 
ents.  To  do  so,  make  use  of  solving  differential  equations  with  the  help  of 
the  Laplace  transformation.  The  transformed  equation  itself  acquires  for  peri 
odic  functions  of  i  the  significance  of  a  frequency  characteristic  when  the 
complex  variable  s  is  substituted  for  by  the  cyclic  frequency  ju>.  However, 
the  transfer  characteristic  is  not  obtained  so  easily  when  periodic  functions 
of  ©  are  assumed  and  the  current  remains  constant.  It  is  therefore  preferable 
to  derive  the  transfer  function  in  the  exact  manner.  The  Laplace  transforma¬ 
tion  of  equation  (4)  for  the  case  ©  =  constant,  then,  is 


U(s)(l+sT)  -  TU^^  =  -MRi(s)  -  MR©si(s)  +  MR©lj.^  .  (15) 

The  initial  conditions  axe  and  i-t_o  =  0-  The  cuxrent  as  a  periodic  func¬ 

tion  :uay  be  expressed  as 

i(t)  =  i  exp(jt--t)  (16) 

and  its  transformation  as 

i!.s)  =  i/(s  -  ju»).  (17) 

The  transformed  differential  equation  then  has  the  appearance; 

U(s)  =  -ITli  (1  +  s©)/(l  +  sT)  (s  -  joj). 


(18) 


By  using  tlie  standardizing  relation  (j),  and.  after  decay  of  the  transi¬ 
ents,  the  time  relation  after  inverse  transformation  is  obtained. 

V  «  (1  +  J‘^)/(l  +  exp(jtdt).  (19) 

This  equation  can  nov  be  discussed.  0  »  T  is  the  matched  condition  and 
for  this  IV  1  »  1.  For  any  other  value  of  the  tine  constants  there  is  a  fre¬ 
quency  dependence  of  V.  In  Fig.  6  the  absolute  value  of  V  is  plotted  against 
the  inverse  frequency  (cycle  length)  in  a  double  logarithmic  scale.  The  ampli¬ 
tude  is  expressed  in  decibels  and  the  circuit  is  matched  when  the  conductivity 

has  a  value  of  10"^^  ohm“lm“^. 


For  cycle  lengths  long  compared  to  2trT  or  2'rT9,  I VI  is  unity  and  no 
errors  exist.  For  cycle  lengths  short  compared  to  2-^1!  or  2tT'9,  |Vl  is  con¬ 
stant  and  has  the  value  9/T.  For  cycle  lengths  between  2:rT  and  2JT9  there  is 
a  smooth  transition  from  unity  to  9/T. 

For  the  case  that  i  is  constant  and  9  is  the  variable,  the  transformation 
of  equation  (4)  is  obtained: 

U/R  (1  +  sT)-  T/R.U^._q  a  -Mi/s  -  sMi  9  (s)  +  Mi9.t^o*  (^O) 

with  U-tsO  =  -MiR,  9  (s)  =  9/(s  - 


Again  normalize  this  equation  and  neglect  the  transients.  Transforming 
this  equation  back  into  the  time  domain, 

V  »  1  +  jw9/(l  +  JdT)  exp(jc^).  (21) 

The  absolute  value  of  the  error  term  is  plotted  in  Fig.  7.  For  periodic 
changes  of  9  which  have  a  cycle  length  long  compared  to  2ir9,  IVl  approaches 
unity  and  there  is  no  influence  of  conductivity  changes  on  current  measure¬ 
ments.  For  a  cycle  length  short  centred  to  2ttT,  an  error  term  of  the  abso¬ 
lute  value  9/T  and  with  the  cyclic  frequency  c*>  is  obtained.  This  error  term 
is  superimposed  on  the  constant  current. 

How  much  conductivity  changes  affect  the  current  recording  in  practical 
cases  cannot  be  stated  definitively  without  knowing  the  average  spectrum  of 
the  conductivity  fluctuations.  However,  examination  of  simultaneous  current 
and  conductivity  recordings  does  not  show  any  clear  transfer  of  conductivity 
fluctuations  on  the  current  reading.  This  is  an  indication  that  the  spectrum 
of  the  conductivity  has  no  significant  component  in  the  cycle  length  region 
short  compeired  to  2n'9. 
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